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Aims. We evaluate the generation of magnetosonic waves in differentially rotating magnetized plasma. 

Methods. Differential rotation leads to an increase of the azimuthal field by winding up the poloidal field lines into the toroidal field lines. An 
amplification of weak seed perturbations is considered in this time-dependent background state. 

Results. It is shown that seed perturbations can be amplified by several orders of magnitude in a differentially rotating flow. The only 
necessary condition for this amplification is the presence of a non-vanishing component of the magnetic field in the direction of the angular 
velocity gradient. 
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\ 1. Introduction 

In astrophysical bodies, differential rotation is often associated 
with magnetic fields of various strength and geometry. If the 
poloidal field has a component parallel to the gradient of angu- 
lar velocity, then differential rotation can stretch toroidal field 
lines from the poloidal ones. In the presence of the magnetic 
field, differential rotation can be the reason for various mag- 
netohydrodynamic instabilities, particularly if the field geome- 
try is complex. Some of these instabilities occur in the incom- 
pressible limit (Velikhov 1959; Fricke 1969; Acheson 1978; 
Balbus & Hawley 1991) which applies if the magnetic field is 
weak and the Alfven velocity is smaller than the sound speed. 
Other instabilities become important only in sufficiently strong 
magnetic fields when the effect of compressibility plays a sig- 
nificant role (Pessah & Psaltis 2005; Bonanno & Urpin 2006, 
2007). Note that incompressible instabilities can often be sup- 
pressed by a strong magnetic field. For example, the magne- 
torotational instability (MRI) does not occur if the magnetic 
field satisfies the condition B 2 > —&TipsQ,Q,'L 2 where p is 
the density, L is the lengthscale of disturbances, 51 = f2(s) 
is the angular velocity, and s is the cylindrical radius; £7' = 
dtt/ds (see, e.g., Urpin 1996, Kitchatinov & Rudiger 1997). 
Moreover, if increases with the cylindrical radius MRI can- 
not arise. 
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In recent years, many simulations of differentially rotating 
magnetized bodies have been performed, and much of the dy- 
namics was interpreted as being a direct consequence of the 
MRI (Brandenburg et al. 1995; Hawley at al. 1995; Matsumoto 
& Tajima 1995; Hawley 2000). Obviously, the MRI cannot 
be the only instability that operates in a rotating magnetized 
gas. For example, stratification can lead to a number of strong 
non-axisymmetric instabilities (Agol et al. 2001; Narayan et al. 
2002; Keppens et al. 2002). Blokland et al. (2005) consider the 
influence of a toroidal field on the growth rate of the MRI and 
find that it leads to overstability (complex eigenvalue). Van der 
Swaluw et al. (2005) study the interplay between different in- 
stabilities and argue that the growth rate of convection can be 
essentially increased due to magnetorotational effects. Note, 
however, that these studies treat the stability of the magnetic 
field with a vanishing radial component, a condition which is 
often not met in astrophysical bodies. In fact, the presence of 
a radial magnetic field can change substantially the stability 
properties (Bonanno & Urpin 2006, 2007). 

In this paper, we consider stability of a differentially rotat- 
ing gas in the presence of a non-vanishing radial magnetic field. 
Differential rotation causes the azimuthal field to increase with 
time by winding up the poloidal field lines into the toroidal 
ones. Therefore, a development of small perturbations occurs 
in the time-dependent background state. We show that stretch- 
ing of the azimuthal field leads to the generation of magne- 
tosonic waves in a flow. Magnetohydrodynamic waves and tur- 
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bulence generated by this instability can play an important role 
in enhancing transport processes in various astrophysical bod- 
ies, such as accretion and protoplanetary Disks, galaxies, stellar 
radiative zones, etc. 

2. Basic equations 

The equations of compressible MHD read 
Vp 1 

v + (v-V)v = — (VxB)xB, (1) 

p Anp 

p + V-(pv) = 0, (2) 

B - V x (v x B) = 0, (3) 

V • B = 0. (4) 

Our notation is as follows: p and v are the gas density and 
velocity, respectively; p is the pressure; B is the magnetic field 
and g is gravity. In this paper, we consider an isothermal gas 
and assume 



p = c p, 



(5) 



where the sound speed, c 2 , is constant. 

We work in cylindrical coordinates (s. (p, z) with the unit 
vectors (e s , e v , e z ). The basic state on which the stability anal- 
ysis is performed is assumed to be axisymmetric with the an- 
gular velocity — Q(s) and B 0. In the presence of the 
non-vanishing radial field B s and differential rotation, the az- 
imuthal field increases with time by winding up the radial field 
lines. If the magnetic Reynolds number is large, then one ob- 
tains from Eq. (4) that the azimuthal field grows linearly with 
time in the basic state, 



B v (t) = B v (0) + sn'B s t, 



(6) 



where B v (0) is the azimuthal field at t = 0. A growth of B v 
given by Eq. (6) can last only while diffusion of the toroidal 
field is negligible. Eventually, in the presence of a finite diffu- 
sivity, a steady state will emerge where winding up is balanced 
by diffusion of the azimuthal field. The time-scale to reach this 
steady state is approximately equal to the diffusion time-scale, 
~ s 2 jr\, where 77 is the magnetic diffusivity. Therefore, our 
consideration is valid only for t < s 2 /rj. However, for astro- 
physical plasmas, this time can be very long because of low 77 
and large length-scale. 

For the sake of simplicity, we assume that gravity is radial, 
g(s) = —g(s)e s , and the basic state is in hydrostatic equilib- 
rium, then 



V« 1 

— = D + — (V xB)xB , 

p 47T/9 



D = g + CI s. 



(7) 



We consider the basic state in which B s oc B v (0) oc B v (t) oc 
1/s. It follows from Eq. (6) that this condition is satisfied if 
sfi' =const. Then, V x B{t) = and the Lorentz force is 
vanishing in Eq. (7). If, additionally, the centrifugal force is 
balanced by gravity, g(s) = s£l 2 , then the pressure and hence, 
density are homogeneous in the basic state. 



Consider stability of axisymmetric short wavelength per- 
turbations with the spatial dependence cx exp(— ik ■ r) where 
k is the wavevector. Small perturbations will be indicated by 
subscript 1, while unperturbed quantities will have no sub- 
script. For the purpose of illustration, we treat in this paper 
the simplest case when k has only the vertical component, 
k = (0, 0, k). Then, the condition V • B 1 = yields B lz = 0. 
The other non-trivial linearized MHD-equations to lowest or- 
der in Ifc • r| _1 read 



dv 



lz 



ikp 



dt 



- + 7— k(Bi ■ B) 
p 4irp 



d P 1 ■ 1 n 
— - zpkv lz = 0, 



dBu 
dt 

dB\ v 
dt 



= iB s kv lz , 

= sfl'Bi s + ikB v (t)v lz . 



(8) 



(9) 



(10) 



(11) 



Since gas is assumed to be isothermal, we have pi/p = pi/p- 
Eliminating all variables from Eqs. (8)-(l 1) in favour of the ver- 
tical velocity perturbation v iz one obtains the following fourth- 
order ordinary differential equation in v\ z : 



dt 4 



where 



k 2 [c 2 s +c\(t)} 



d 2 v lz 



dt 2 

'„2 2 



6w 3 



BiY 



dvu 
dt 



+6k"c" As (sn') 2 v lz =0, 



(12) 



CA v (t) = B v {t)/^/AKf> , u>% n = k 2 c A(p (t)c As stt'. 

Eq. (12) with the corresponding initial conditions describes 
the behaviour of velocity perturbations on the time-dependent 
background state. The behaviour of perturbations of the ra- 
dial and azimuthal magnetic field can then be calculated from 
Eq. (10) and (11). 



3. Numerical results 

To follow the behaviour of perturbations, it is convenient to 
introduce the dimensionless quantities 



t = Qt , q 



kH , H 



CA<pO 



B V (P) 



Then, introducing v = v\ z (t) jv\ z (0) where v\ z (0) is the initial 
velocity perturbation, we obtain from Eq. (12) 

d 4 v 9r , . , r—- 1—- ,o,d 2 v 



dv 



6qafy/0 a (y/0 v + qy* 0,t)— + 6q 2 x 2 f3 s v = 

(IT 



(13) 



The dependence of the solution on the wavelength is character- 
ized by the parameter x and, on the initial magnetization of gas, 
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Fig. 1. The time dependence of v (solid line), C'a s (dashed 
line), and Ca v (dashed-and-dotted line) for q = 0.1, x = 3, 
(i s = 0.01, and /3 v = l. 



by the parameters (3 S and (3 V . To solve Eq. (13), one needs the 
initial conditions for three time derivatives of v. In actual cal- 
culations we try different initial conditions because their choice 
is determined by the origin of perturbations which is uncertain. 

Eqs. (10)-(1 1) can also be written in a dimensionless form. 
We have 



dCAs 

dr 



ixV/3 s w(r), 



dCA v 

dr '' 
where 

Cas = 



qC As + ix{->/frp + q\/%T)v(T), 



B u {t) 



Cau 



Bi v {t) 



(14) 
(15) 

(16) 



In calculations, we assume that the initial perturbations of the 
magnetic field are vanishing. Eqs. ( 1 3)-( 15) were solved nu- 
merically for a wide range of the parameters. 

In Fig. 1, we plot the time dependence of v, Ca s , and Ca v 
in the regions where the angular velocity increases slowly with 
s (q = 0.1). The other parameters are x — 3, (i s = 0.01, and 
V = 1. Integrating Eq. (13), we choose the following initial 
conditions: dv(0)/dr = 1 and d 2 v{0)/dT 2 = d 3 v(0)/dr 3 = 
0. The behaviour of the velocity and magnetic field turns out to 
be substantially different. During the initial stage, the velocity 
reaches a flat maximum at r <~ 60 — 70 which corresponds 
approximately to 10P where P is the rotation period. At that 
time, the perturbation of velocity is w 25 times greater than 
its initial value. Perturbations of the magnetic field are much 
larger. After r ~ 100 (t ~ 15P), Ca s approaches the satu- 
ration level that is approximately a factor 10 3 greater than the 



Fig. 2. Same as in Fig. 1 but for q = —0.1 and x = 10. 



initial perturbation of v\ z . The amplification of the azimuthal 
field is even larger. After t ~ 15P, Ca v is a factor 10 4 stronger 
than the initial perturbation vi z (0), and it continues to grow 
approximately linearly with time. Note that the long-term be- 
haviour (r 3> 200) of these quantities is also different. The ve- 
locity v changes the sign and exhibits some sort of oscillatory 
behaviour with the amplitude growing oc t 1 / 2 whereas Ca s 
reaches the saturation level. A perturbation of the azimuthal 
field Bi v is still oc t at very large t. 

In contrast with the magnetorotational instability, the be- 
haviour of perturbations is not crucially sensitive to the sign 
of O'. To illustrate this point, we show in Fig. 2 the time de- 
pendence in the case q = —0.1 (51 decreases slowly with s). 
Qualitatively, the dependences are the same. The only essen- 
tial difference is that the azimuthal field generated by differen- 
tial rotation has a negative sign because the radial gradient of 
the angular velocity is negative. Initially (r < 100), the main 
contribution to B\ v is provided by compressibility (the second 
term on the r.h.s. of Eq. (11)). Later on, however the azimuthal 
magnetic perturbations are determined by winding up the ra- 
dial field perturbations and becomes negative. In Fig. 2, the 
modulus of Ca v is shown. Since x is larger in this figure, per- 
turbations of all quantities can reach higher values compared to 
Fig. 1. For instance, the velocity becomes ~ 200 times greater 
than the initial value after 15 rotation periods. The radial pertur- 
bation of the Alfven velocity reaches the saturation level which 
is <~ 2 x 10 4 times greater than v\ z (0), and the azimuthal pertur- 
bation is even about 10 times larger. For very long t(r » 200), 
the velocity exhibits the oscillatory behaviour with the ampli- 
tude oc t 1 / 2 . In Fig. 3, we plot the behaviour of perturbations 
in the case of a very weak differential rotation, q = 0.01. For 
a weak rotation shear, the amplification of initial perturbations 
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Fig. 4. Same as in Fig. 1 but for q = 1 and f3 s = (3 V = 0.1. 



is slower. The velocity reaches the maximum after <~ 40P, and 
this maximum is rather high: v\ z exceeds its initial value by a 
factor ~ 100. Despite a slower growth rate, saturation of Ca s 
still occurs at a rather high level (Ca s ~ 2x 10 4 ). As expected, 
the azimuthal field grows linearly with time and can reach a 
very high value. 



The same dependences, for a relatively large value of q, 
are shown in Fig. 4. In this case, perturbations exhibit a clear 
oscillatory behaviour. This concerns particularly perturbations 
of the vertical velocity. After a short initial relaxation (last- 
ing ^2 — 3 rotation period), v reaches a self-similar regime 
when the amplitude of oscillations grows with time oc r 1 / 2 
and the frequency oc r. Perturbations of the magnetic field also 
exhibit a weak oscillatory behaviour but the amplitude of oscil- 
lations is much smaller. This is because the field components 
are given by the integrals of a rapidly oscillating function v 
(see Eq. (10) and (11)), and oscillations in E>i s and B\ v are 
smoothed. Perturbations of the radial field reaches a saturation 
level after <~ 2 — 3P. In saturation, the radial Alfven velocity 
Cas is approximately 10 times greater than the initial pertur- 
bation of the vertical velocity. Perturbations of the azimuthal 
field grow rather rapidly at the initial stage however, later on 
(at t > 2P), they continue to grow linearly with time. 

In Fig. 5, we show the long-term behaviour of the ampli- 
tude of the vertical velocity. Except for the initial stage, the nu- 
merical results are well fitted by the same dependence oc r 1 / 2 
as that in Fig. 4. We can compare these numerical results with 
the analytical solution of Eq. (13) for large r. It can be easily 
checked that the asymptotic solution of Eq. (13) at large r is 
given by 



v oc Vr cos (^qx^ffT s T 2 + ^ + 0(t~ 



3/2 



)■ 



(17) 



Indeed, it can be noted that the vertical velocity is oc r 1 / 2 
and the frequency of oscillations is oc r at r — > oo in com- 
plete agreement with our numerical results. The frequency of 
oscillations depends on the magnetic field and is higher for a 
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stronger initial radial field. Note that the asymptotic solution 
(17) is valid for a wide range of parameters and initial condi- 
tions. It is clear from solution (17) that Bi s should reach satu- 
ration quite rapidly. We have from Eq. (14) that 

C As (t) = ix^/Js I v(r')dT'. (18) 
Jo 

The time integral of the asymptotic part of v (Eq. (17)) gives 
a negligible contribution to Ca s , and the only non-vanishing 
contribution can be provided by the value of v at relatively 
short r. If Cas reaches saturation rapidly then, as it follows 
from Eq. (15), Ca v should grow approximately a rat large r. 
This sort of behaviour agrees perfectly with the results of our 
numerical calculations. 



time integrals of a rapidly oscillating velocity. In the case of 
a strong differential rotation (q <~ 1), perturbations exhibit the 
oscillatory behaviour from the very beginning and the initial 
growth of their amplitude is less significant. 

The generation of magnetosonic waves occurs even if the 
magnetic field is very strong and suppresses different MHD- 
instabilities which can arise in a differentially rotating flow (for 
example, the MRI). The presence of differential rotation and 
radial magnetic field is, however, crucially important for the 
considered process. Since both differential rotation and radial 
field are quite common in astrophysics, we believe that the con- 
sidered mechanism can occur in various astrophysical bodies 
and plays an important role in enhancing transport processes in 
plasma. 



4. Discussion 

We have shown that the winding up of toroidal field lines from 
the poloidal field lines is accompanied by an amplification of 
seed perturbations of the velocity and magnetic field. A very 
simplified model has been considered in this paper, but we 
believe that qualitatively the same results can be obtained for 
more general background states and perturbations. Differential 
rotation and compressibility of the gas lead to the generation of 
magnetosonic waves with the amplitude that grows with time. 
The physical processes responsible for this amplification are 
exactly the same that result in the instability considered by 
Bonanno & Urpin (2006, 2007). The only difference is that, 
in this paper, we consider the development of perturbations on 
a time-dependent background state and, as a result, the growth 
of perturbations is not exponential. 

The behaviour of seed perturbations depends essentially 
on various parameters and can generally be rather complex. 
If the parameter q — sfi'/fi is relatively small (< 0.1), then 
the perturbations of velocity and magnetic field initially grow 
monotonously and can reach quite high values. For example, 
the perturbation of the vertical velocity becomes approximately 
150 — 200 times greater than its initial value after only 15-30 
rotation periods (see Figs. 2 and 3). Perturbations of the mag- 
netic field reach even higher values during the initial stage. For 
instance, the Alfven velocity corresponding to the perturbation 
of the radial field component, Ca s , can exceed the initial veloc- 
ity perturbation by a factor ~ (1 — 3) x 10 5 after the same time, 
but the perturbations of the toroidal field are even stronger. As 
a result of such a strong initial amplification, seed perturbations 
can already reach a non-linear regime after 15-30 rotation peri- 
ods if their initial values are sufficiently large. Further evolution 
of perturbations will then be entirely determined by non-linear 
effects. However, if the non-linear regime is not reached during 
this initial stage, the behaviour of perturbations becomes oscil- 
latory with slowly growing amplitude (oc i 1 / 2 ). At sufficiently 
large t, the frequency of oscillations grows linearly with time 
and is given approximately by 

lj w n (^kHy/fasn'tj . (19) 

Note that perturbations of the magnetic field exhibit more reg- 
ular behaviour because they can be expressed in terms of the 



Acknowledgments. This research project has been supported 
by a Marie Curie Transfer of Knowledge Fellowship of the 
European Community's Sixth Framework Programme un- 
der contract number MTKD-CT-002995. VU thanks INAF- 
Ossevatorio Astrofisico di Catania for hospitality. 

References 

Acheson D. 1978. Phil. Trans. R. Soc. London, 289A, 459 

Agol E., Krolik J., Turner N., Stone J. 2001. ApJ, 558, 543 

Balbus S.A., & Hawley J.F. 1991. ApJ, 376, 214 

Blokland J.W.S., van der Swaluw E., Keppens R., Goedbloed J.P. 

2005. A&A, 444, 337 
Bonanno, A. & Urpin, V. 2006. Phys. Rev. E, 73, 066301 
Bonanno, A. & Urpin, V. 2007. ApJ, 662, 851 

Brandenburg, A., Nordlund, A., Stein, R., Torkelsson, U, 1995, ApJ, 

446, 741 
Hawley, J. 2000. ApJ, 528, 462 

Hawley, J., Gammie, C, & Balbus, S. 1995. ApJ, 440, 742 
Fricke K. 1969. A&A, 1, 388 

Keppens, R., Casse, E, & Goedbloed, J. 2002, ApJ, 569, L121 
Kitchatinov, L.L., & Rudiger, G. 1997, MNRAS, 286, 757 
Matsumoto, R., & Tajima, T. 1995, ApJ, 445, 767 
Narayan R., Quataert E., Igumenshchev I., Abramowicz M. 2002. 

ApJ, 577, 295 
Pessah, M. & Psaltis, D. 2005. ApJ, 628, 879 
Urpin, V. 1996. MNRAS, 280, 149 

Van der Swaluw E., Blokland J.W.S., & Keppens R. 2005. A&A, 444, 
347 

Velikhov, E.P. 1959. Sov. Phys. JETP, 9, 995 



